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Abstract 

Using a generalization of vector calculus for the case of non-integer dimensional space 
we consider a Poiseuille flow of an incompressible viscous fractal fluid in the pipe. Fractal 
fluid is described as a continuum in non-integer dimensional space. A generalization of 
the Navier-Stokes equations for non-integer dimensional space, its solution for steady flow 
of fractal fluid in a pipe and corresponding fractal fluid discharge are suggested. 


1 Introduction 

A cornerstone of fractal fluids is the non-integer dimension EM- The mass of fractal fluid 
satisfies a power law relation M ~ R D , where M is the mass of the ball region with radius R , 
and D is the mass dimension [4], Fractal fluid can be described by four different approaches: 
(1) Using the methods of’’Analysis on fractals” [5| [6], U, 8];.9J[X0] it is possible to describe fractal 
media; (2) An application of fractional-differential continuum models suggested in [TTl 12] , and 
then developed in [13], 111] [15l 16] [17] E] , where so-called local fractional derivatives [19] are 
used; (3) Applying fractional-integral continuum models suggested in [20j |2H EZ, 23] [4] (see 
also |24|-[28| and ]2§]-|32j). where integrations of non-integer orders and a notion of density of 
states [4] are used; (4) Fractal media can be described by using the theory of integration and 
differentiation for a non-integer dimensional space [33, Si] [35] (see also mm- 

Let us note that main difference of the continuum models with non-integer dimensional 
spaces form the fractional continuum models suggested in [20] 2T] [22] 23] 0J may be reduced 
to the following, (a) Arbitrariness in the choice of the numerical factor in the density of 
states is fixed by the equation of the volume of non-integer dimensional ball region, (b) In 
the fractional continuum models suggested in [20] 21,1], the differentiations are integer orders 
whereas the integrations are non-integer orders. In the continuum models with non-integer 
dimensional spaces the integrations and differentiations are defined for the spaces with non¬ 
integer dimensions. 

In this paper, we consider approach based on the non-integer dimensional space. The power 
law M ~ R d can be naturally derived by using the integrations in non-integer dimensional 
space [33], where the mass dimension of fractal fluid is connected with the dimension of this 
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space. A vector calculus for non-integer dimensional space proposed in this paper allows us 
to use continuum models with non-integer dimensional spaces to describe for fractal fluids. 
This is due to the fact that although the non-integer dimension does not reflect completely 
the geometric properties of the fractal media, it nevertheless permits a number of important 
conclusions about the behavior of fractal structures. Therefore continuum models with non¬ 
integer dimensional spaces can be successfully used to describe fractal fluids. 

Integration over non-integer dimensional spaces are actively used in the theory of criti¬ 
cal phenomena and phase transitions in statistical physics [36l 37] , and in the dimensional 
regularization of ultraviolet divergences in quantum field theory [ 38|[39l 53]. The axioms for in¬ 
tegrations in non-integer dimensional space are proposed in [40J [34J and this type of integration 
is considered in the book by Collins [33] for rotationally covariant functions. In the paper [34] a 
mathematical basis of integration on non-integer dimensional space is given. Stillinger [33] sug¬ 
gested a generalization of the Laplace operator for non-integer dimensional spaces also. Using 
a product measure approach, the Stillinger’s methods [33] has been generalized by Palmer and 
Stavrinou [35] for multiple variables case with different degrees of confinement in orthogonal 
directions. The scalar Laplace operators suggested by Stillinger in [34] and Palmer, Stavrinou 
in [33] for non-integer dimensional spaces, have successfully been used for effective descriptions 
in physics and mechanics. The Stillinger’s form of Laplacian for the Schrodinger equation in 
non-integer dimensional space is used by He [4T] 02J 43] to describe a measure of the anisotropy 
and confinement by the effective non-integer dimensions. Quantum mechanical models with 
non-integer (fractional) dimensional space have been discussed in [34, 33] 04, [43] 4B] 07], 43] 
and gSHSHEUGS]. Recent progress in non-integer dimensional space approach also includes 
description of the fractional diffusion processes in non-integer dimensional space in [[53], and the 
electromagnetic fields in non-integer dimensional space in [34, [331 ESJ and [57] [58, [59] [60, 61]. 

Unfortunately, [34] [33] proposed only the second order differential operators for scalar fields 
in the form of the scalar Laplacian in the non-integer dimensional space. A generalization of 
the vector Laplacian p2J for the non-integer dimensional space is not suggested in [34] [35 ]. 
The first order operators such as gradient, divergence, curl operators are not considered in 
[341 [33] also. In the work [61J the gradient, divergence, and curl operators are suggested only as 
approximations of the square of the Laplace operator. Consideration only the scalar Laplacian 
in the non-integer dimensional space approach greatly restricts us in application of continuum 
models with non-integer dimensional spaces for fractal fluids and material. For example, we 
cannot use the Stillinger’s form of Laplacian for vector field v(r, £) in hydrodynamics of fractal 
fluids, in fractal theory of elasticity and thermoelasticity, in electromagnetic theory of fractal 
media to describe processes in the framework non-integer dimensional space approach. 

In this paper, we propose to use a vector calculus for non-integer dimensional space, and we 
define the first and second orders differential vector operations such as gradient, divergence, the 
scalar and vector Laplace operators for non-integer dimensional space. In order to derive the 
vector differential operators in non-integer dimensional space we use the method of analytic 
continuation in dimension. For simplification we consider rotationally covariant scalar and 
vector functions that are independent of angles. It allows us to reduce differential equations in 
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non-integer dimensional space to ordinary differential equations with respect to r. The proposed 
operators allow us to describe fractal media to describe processes in the framework of continuum 
models with non-integer dimensional spaces. In this paper we describe a Poiscuillc flow of an 
incompressible viscous fractal fluid in the pipe. A generalization of the Navier-Stokes equation 
for non-integer dimensional space to describe for fractal fluid are suggested. A solution of this 
equation for steady flow of fractal fluid in a pipe and corresponding fractal fluid discharge are 
derived. 


2 Fractal fluids 

A basic characteristic of fractal fluids is the non-integer dimensions such as mass or ’’particle” 
dimensions [3]. For fractal fluids the number of particles Nd(W ) or mass Md(W) in any region 
W C M 3 of this fluid increase more slowly than the 3-dimensional volume V 3 (W) of this region. 
For the ball region W with radius R in an isotropic fractal fluid, this property can be described 
by the relation between the number of particles Nd(W ) in the region W of fractal fluid, and 
the radius R in the form 


N d (W) = N 0 (R/R 0 ) d , R/R 0 » 1, (1) 

where i? 0 is the characteristic size of fractal fluid such as a minimal scale of self-similarity of 
a considered fractal fluid. The number D is called the ’’particle” dimension. It is a measure 
of how the fluid particles fill the space. The parameter D does not depend on the shape of 
the region W. Therefore fractal fluids can be considered as fluid with non-integer ’’particle” or 
mass dimension. 

If the fractal fluid consists of particles with identical masses mo, then relation ([!]) gives 

M d (W) = M 0 (R/R 0 ) d , R/R 0 » 1, (2) 

where M 0 = m 0 N 0 . In this case, the mass dimension coincides with the ’’particle” dimension. 

As the basic mathematical tool for continuum models of fractal fluids, we propose to use the 
integration and differentiation in non-integer dimensional spaces. In Section 7, we will show 
that the power-law relation (J2]) for an isotropic fractal fluid can be naturally derived by using 
the integration over non-integer dimensional space, where the space dimension is equal to the 
mass dimension of fractal fluid. 

In order to describe fractal fluid by continuum models with non-integer dimensional spaces, 
we use the concepts of density of states c 3 (D,r) that describes how closely packed permitted 
places (states) in the space M 3 , where the fractal fluid is distributed. The expression dVoijc) = 
c 3 (D,r)dV 3 is equal to the number of permitted places (states) between V 3 and V 3 + dV 3 in 
M 3 . The notation d D r also will be used instead of dVr, >(r). Note that density of states and 
distribution function are different concepts, and it is impossible to describe all properties of 
fractal fluids by the distribution function only. 
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For fractal fluids, we can use the equation 

dN D {W) =n(r)dVb(r), (3) 

where n(r) is a concentration of particles that describes a distribution of number of particles on 
a set of permitted places (possible states). The density of states is chosen such that dVo( r) = 
Cs(D, r) dV 3 describes the number of permitted states in dV 3 . 

The form of the function c${D,y) is defined by symmetries of considered problem and 
properties of the described fractal fluid. A general property of density of states for fractal 
fluids is a power-law type of these functions that reflects a scaling property (fractality) of the 
fractal fluid. To simplify our consideration in this paper we will consider only isotropic fractal 
fluids with density of states that is independent of angles. In this case, the form of density of 
states is defined such that dVo is an elementary volume of the non-integer dimensional space. 

In the continuum models of fractal fluids, we should work with the dimensionless variables 
x/Rq —> x, y/Ro —> x, z/R 0 —> x, r /R 0 — > r, in order to physical quantities of fractal fluids 
have correct physical dimensions. 

3 Vector differential operators in non-integer dimensional 
space 

To derive equations for vector differential operators in non-integer dimensional space, we use 
equations for the differential operators in the spherical (and cylindrical) coordinates in R n for 
arbitrary n to highlight the explicit relations with dimension n. Then the vector differential 
operators for non-integer dimension D can be defined by continuation in dimension from integer 
n to non-integer D. To simplify we will consider only scalar fields (p and vector fields v that 
are independent of angles 

<p(r) = <p(r), v(r) = v(r) = v r e r , 

where r = |r| is the radial distance, e r = r/r is the local orthogonal unit vector in the directions 
of increasing r, and v r = v r (r) is the radial component of v. We will work with rotationally 
covariant functions only. This simplification is analogous to the simplification for definition of 
integration over non-integer dimensional space described in Section 4 of the book [33] . 

3.1 Vector differential operators for spherical and cylindrical cases 

Using the continuation from integer n to arbitrary non-integer D, we can get explicit definitions 
of differential operators for non-integer dimensional space in the following forms. Note that 
the same expressions can be obtained by using the integration in non-integer dimensional space 
and the correspondent Gauss’s theorem [72] 

Let us define the differential vector operations such as gradient, divergence, the scalar and 
vector Laplacian for non-integer dimensional space. For simplifications, we assume that the 
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vector field v = v(r) be radially directed and the scalar and vector fields <^(r), v(r) are not 
dependent on the angles. 

The divergence in non-integer dimensional space for the vector field v = v(r) is 


Div?v = ar + 


dv r D 


V r - 


The gradient in non-integer dimensional space for the scalar held p = p{r) is 


Grad? p = ^ e r . 
or 


(4) 


(5) 


The scalar Laplacian in non-integer dimensional space for the scalar held p = p(r) is 


SAjy=Div?Grad ??=*? + £=!* 


( 6 ) 


The vector Laplacian in non-integer dimensional space for the vector held v = v(r) e r is 

'd 2 v r D — 1 dv r D — 1 


v a d 


A?v = Grad? Div? v = 


dr 2 


+ 


dr 


„2 Vr I er " 


(7) 


If D = n, equations (jU17]) give the well-known formulas for integer dimensional space M n . 
Let us consider a case of axial symmetry of the huid, where the helds <p(r) and v(r) = v r (r) e r 
are also axially symmetric. We will direct the Z-axis along the axis of symmetry. Therefore we 
use a cylindrical coordinate system. 

The divergence in non-integer dimensional space for the vector held v = v(r) is 


Div? v = — +- v r . 

or r 

The gradient in non-integer dimensional space for the scalar held ip = <p(r) is 

n a D ^ 

Grad p = — e r . 

or 

The scalar Laplacian in non-integer dimensional space for the scalar held p = p(r) is 

° 2 p D- 2 dp 


dv r D — 2 


( 8 ) 


(9) 


’A ?p = 


2 Qrp 


( 10 ) 


The vector Laplacian in non-integer dimensional space for the vector held v = v(r) e r is 


v\D 


A? v = 


d 2 v r D — 2 dv r D — 2 

H---7T— v r ) e r . 


dr 2 


r dr 


( 11 ) 
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Equations ffSl fTO) can be easy generalized for the case p = p(r, z ) and v(r, z) = v r (r, z) e r + 
v r (r, z) e z . In this case the curl operator for v(r, z) is different from zero, and 


Curlf v 


dv r dv z \ 

dz dr ) e 


( 12 ) 


For D — 3 equations O - CD and (jT 2 T) give the well-known expressions for the gradient, 
divergence, curl operator, scalar and vector Laplacian operators 

The suggested operators for 0 < D < 3 allow us to reduce .D-dimensional vector differen¬ 
tiations (HD - (HD and © - (HID to derivatives with respect to r = |r|. It allows us to reduce 
partial differential equations for fields in non-integer dimensional space to ordinary differential 
equations with respect to r. 


3.2 Stillinger’s Laplacian for non-integer dimensional space 

For a function p — p(r, 6) of radial distance r and related angle 0 measured relative to an axis 
passing through the origin, the scalar Laplacian in a non-integer dimensional space proposed 
by Stillinger [33] is 


St^D 


1 d ( D _ 1 d \ 1 d ( . D _2 n d \ 

r n -' dr (' fr) + r 2 sin D - 2 »S» (™ $ &>) ’ 


(13) 


where D is the dimension of space (0 < D < 3), and the variables r > 0, 0 < 6 < n. Note that 
( Si A D ) 2 ^ st A 2D . If the function depends on radial distance r only (p = p(r)), then 


St A D p(r) 


1 9 Ld- i dvirY \ 
r D ~i Q r dr J 


d 2 p(r) + D — 1 dp(r) 


(14) 


It is easy to see that the Stillinger’s form of Laplacian st A D for radial scalar functions <p(r) = 
<p(r) coincides with the scalar Laplacian S A® defined by (ED ? i- e -> 

st A D ip(r) = s A D tp(r). (15) 


The Stillinger’s Laplacian can be applied for scalar fields only. It cannot be used to describe 
vector fields v = v r (r) e r because this Laplacian for D = 3 is not equal to the usual vector 
Laplacian for R 3 , 

. ■, . . , . . . (d 2 v r 2 dr v 2 \ . . 

A v(r) + A v(r) = + ~ ~ ^ v r ) e r . (16) 

The gradient, divergence, curl operator and vector Laplacian are not considered by Stillinger 
in paper [34] . 
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3.3 Differential operators for d ^ D — 1 

Let us consider a ball region B D in the fractal fluids with the boundary Sj, = dB D with 
dimensions 

dim(U D ) = D , dim(5d) = d. (17) 

Equations (j3J - (TTIT) define differential operators for spaces with non-integer dimension D and 
boundary dimensions d — D — 1. In general, the dimension D of the region of a fractal fluid 
and the dimension d of boundary of this region are not related by the equation d — D — 1, i.e., 

dim (dB D ) ^D- 1. (18) 


Using the integration in non-integer dimensional space and the correspondent Gauss’s the¬ 
orem, we can define the divergence for the case d ^ D — 1 by the equation 


Div 


DA 


V = 


^d+i-D)^ y(D/2) 


T((d+l)/2) \r 


oD—l—d 


dvJr) d . . 


(19) 


For d — D — 1, we get (Ul) . We can define the parameter 


a r = D — d, (20) 

that can be interpreted as a dimension of fractal fluid along the radial direction e r . Using (l20j) . 
equation (1191) can be rewritten in the form 


Div?’ d v = 7 j-( 1—Qr ) /2 


r((d + cn r ) / 2) 
r((d+l)/ 2 ) 


1 dv r (r ) 

rjrtOLj' 1 



( 21 ) 


This is divergence operator for non-integer boundary dimension d ^ D — 1. For a r = 1, we 
have d = D — 1, and then equations (fT9|) . f[2Tj) give (J4|). 

The gradient for the scalar field <p(r) = <p(r) and the radial dimension a r ^ 1 is defined by 


Grad ?> d <p 


r(cx r /2) dip{r) 

j2 rj^OL'p 1 ^) r f‘ 


( 22 ) 


For oi r = 1, equation (1221) gives (T3T) . 

Using the operators (|22l) and ffl9|) for the fields p> = <p(r) and v = v(r) e r , we can get the 
scalar and vector Laplace operators for the case d ^ D — 1 by 


5 Af’V = Div^Grad^’^, 


F A^v = Grad?’ d Div^v. 


(23) 


The scalar Laplacian for d ^ D — 1 for the field ip = <p(r) is 


s A?’ d cp = A(d,a r ) 


1 d 2 ip d + 1 — a r d(p 


r 2a r 2 


dr 2 




dr 


(24) 
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where 


(25) 


A(d, a r ) = 


r [yd + 0L r )/ 2) r(o' r /2) 


7 r “r-l/2 r ((d+ i)/2) ' 
The vector Laplacian for d 7 ^ .D — 1 for the field v = v(r) e r is 


v A?’ d v = A(d, Or) f - 


1 d 2 v r d + 1 — a r dv r d 


■2a r -2 


dr 2 


+ 


a r 


«2oif’ X 


dr 


r’2ctj’ 


V r I 


(26) 


The vector differential operators (j22J), (TTTJp . (12X1) and (1251) allow us to describe complex 
fractal fluids with the boundary dimensions d 7 ^ D — 1 by continuum models with non-integer 
dimensional spaces. 


4 Navier-Stokes equations in non-integer dimensional space 
for fractal fluid 

A motion of an incompressible viscous fractal fluid in the framework of continuum model with 
non-integer dimensional space is described by the equations 


Div ) 9 v = 0, 

— f — - Grad ) 9 p + v 1 A ) 9 v, 
at p 


(27) 

(28) 


where f is the vector held of mass forces, v is the kinematic viscosity is the ratio of the dynamic 
viscosity /1 to the density of the fluid p, and d/dt is the material derivative 


dv <9v 
dt dt 


v, Grad ) 9 v 


(29) 


In equations fl2T]) -([29 ]) the gradient Grad) 9 , the divergence Div) 9 , and the vector Laplacian v A)? 
are defined by equations ©, (j4j), ([7]) for spherical symmetry. For cylindrical symmetry, these 
operators are defined by equations ©, (ed, m- 

If the dimension D of the region of a fractal fluid and the dimension d of boundary of this 
region are not related by the relation d = D — 1, i.e., a r — D — d ^ 1, then we should use the 
equations 


dv 

dt 


= f 


Divf’ d v = 0, 

(30) 

-Gradf’ d p + z/ y A^ d v, 

P 

(31) 

^ + (v, Gradfv), 

(32) 


where the gradient Grad) 9 ^, the divergence Div) 9 ^, and the vector Laplacian v A^ ,d are defined 
by equations (| 22 |) . fl 2 ip . (I26|) . 
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Equations (1251) and (|3T1) can be called the Navier-Stokes equations for non-integer dimen¬ 
sional space. 

It is convenient to work in the dimensionless space variables x/Rq — > x, y/Ro — > x, z/Rq — > 
x, r/Ro —> r, that yields dimensionless integration and dimensionless differentiation in non¬ 
integer dimensional space. Here i ? 0 is the characteristic size of a fractal fluid, which is always 
finite. For example, i ? 0 can be the minimal scale of self-similarity of a considered fractal fluid. 
Then the density is properly scaled such that the mass Q of fractal fluid and the fields v, p, f 
have correct physical dimensions. 

The Navier-Stokes equations (12Hj) and (l3Tj) describe dynamics of fractal fluids in the frame¬ 
work of continuum models with non-integer dimensional spaces. These equations allow us to 
describe the isotropic fractal fluid only when the presence of spherical or cylindrical symmetry. 

Equations (1251) and (ITT]) can be used, when the fields p, v, v have the form p = p(r) and 
v = v r (r ) e r , f) = f r (r ) e r does not depend on the angles. 

As an example of application of the Navier-Stokes equations ([25]) and (l3T|) . we consider 
a steady flow of fractal fluid in a pipe, and fractal fluid discharge in the next sections. In 
the next sections, we derive the Poiseuille equation for fractal fluids from the Navier-Stokes 
equations (125|) and (T5Tj) . Equations can be used for any other problems of hydrostatics and 
hydrodynamics of fractal fluids within the case of spherical and cylindrical symmetries 

To consider anisotropic fractal fluids, and problems without spherical and cylindrical sym¬ 
metries, we cannot use the Navier-Stokes equations (l28|) and (l3Tj) . In this case, we should apply 
a product measure approach [73]. The product measure approach to describe fractal properties 
of space-time has been considered in |63j . The product measure approach for the fractional 
spaces has been suggested in [6H IS], where fractional phase space is considered with its inter¬ 
pretation as a non-integer (fractional) dimensional space. For non-integer dimensional spaces 
the product measure approach is suggested in [35], where each orthogonal coordinates has own 
dimension. The product measure approach for the fractional-integral continuum models has 
been considered in [25, 2BJ127, 281 • 


5 Steady flow of fractal fluid in a pipe 

In this section, we derive the Poiseuille equation for fractal fluids from the Navier-Stokes equa¬ 
tions (125|) . Let us consider a simple problem of motion of an incompressible viscous fractal 
fluid. Using the continuum models with non-integer dimensional space, we describe a steady 
flow of fractal fluid in a pipe with circular cross-section. We take the axis of the pipe as the 
X-axis. The velocity for laminar of fractal fluid is along the A"-axis at all points, and is a 
function of r only 

v = v(r) = v x (r)e x . (33) 

The equation of continuity is satisfied identically. The components of the Navier-Stokes equa¬ 
tion for E-axis and Z -axis give that the pressure is constant over the cross-section of the pipe. 
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We shall solve the equation for a pipe with circular cross-section. Taking the origin at 
the center of the circle we can use cylindrical symmetry v = v(r) = v x (r)e x . Using the 
Navier-Stokes equations fl28l) . we have 


s 


A ?v x (r) 


1 dp 
p dx ’ 


(34) 


where p = pis, and dp/dx is a constant. The pressure gradient dp/dx may be written — Ap/l , 
where A p is the pressure difference between the ends of the pipe and l is its length. 

Using the scalar Laplacian for non-integer dimensional space the Navier-Stokes equation 
()34l) takes the form 


d 2 v x (r) D — 2 dv x 1 dp 

dr 2 r dr p dx 


(35) 


For 1 < D < 3 and 0 < D < 1, the general solution of (1351) is 


v x (r 


= Ci r 3 ~ D + C 2 + 


2 (D — 1) p dx 


dp 2 
r 


(0 < D < 3, D^l). 


(36) 


For D 

For D 


3, 

1 , 


we have 

v x (r) = Ci ln(r) +C 2 + 
we get the general solution 


J_dp_ r 2 

Apdx 


(37) 


v x (r) = Ci r 2 + C 2 + r 2 (2 ln(r) - 1 ) . (38) 

4/i dx 

It should be noted that dimensions D — 1 of the fractal fluid do not correspond to the distri¬ 
bution of particles along the line. The fractal media with D = 1 describe a distribution of fluid 
particles in 3-dimensional space such that the mass dimension of the distribution is equal to 
D = 1. 

Let us determine a flow of fractal fluid in a pipe of annular cross-section with the internal 
radius Ri and external radius R 2 . The constants Ci and C 2 in the general solution are 
determined from the boundary conditions 


? 

'S 

to 

O 

(39) 

Using these conditions have the form 


Ci r\~ d + Co + , 1 , d 7 p Ri = o, 

11 ^ 2(D-l)fj,dx 1 

(40) 

CiR 3 2 ~ d + C 2 + , 1 dp R\ = 0. 

12 2(D — 1) pdx 2 

(41) 
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Then the constants are 


1 


(42) 


Cl = 


dp 


R\ 


R\ 


2 (D — 1) n dx R:'\ 


3-D 


R 


3-D ’ 


Co — 


dp R\~ d R\ - R\~ d R\ 


Substitution of 


and 

1 


v x r = 


2 (D — 1) p dx 
into gives 

dp ( R 2 - R 2 


R\ 


3-D 


R 


3-D 


2 (D — 1) p dx \R\ 


j3-D 


where 0 < D < 1 and 1 < D < 3. 
Using the variables 


x = 


R: 


Ro 

Ri 


3-D 


TOO 

r 3 -D + 


3-D 


R\ 


R: 


3-D 


Ri 


3-D 


R\ 


3-D 


R\ 2 

— + r 2 


r 

y = TTi 


equation dTil) can be represented in the form 


v(x,y) = 


Rl 


dp 


2 {D — 1) p dx \ 1 


1 — x 2 


x ' 


3-D 


Ri 


y 3 ~ D + 


x 2 _ X 3-D 


X 


3-D 


+ y‘ 


(43) 


(44) 


(45) 


To demonstrate some properties of the velocity v x {r ) defined by (1441) . we can visualize the 
function (l45lh for x G [1; 100], y G [1; 100] and different values of dimensions D = 2.9, D = 2.7, 
D = 2.0, D = 1.1. The plots of function are presented by Figures 1-4, where pi — 1, R\ — 1 
and dp/dx = —1. 

The flow in a pipe of annular cross-section with the radius R , i.e. R± = 0 and Ro = R , we 
have 


(r) = 


2 {D — 1) p dx 


d P R 2 


3-D 


r \ 2 

RJ 


(0 < D < 3, D^l). 


(46) 


Equation (146|) can be called the Poiscuillc equation for flow of fractal fluid. For the case of 
non-fractal fluid (D — 3), equation (1461) gives the well-known Poiscuille equation 


v x (r) 


4 pi dx \ 



(47) 


Thus the velocity distribution across the pipe is parabolic for the non-fractal fluids. For the 
fractal fluids, we have non-integer power-law (l46|h 

Note that suggested Poiscuillc equation for fractal fluids, which are derived from the pro¬ 
posed Navier-Stokes equation for non-integer dimensional space can be used only to describe 
flow of fractal fluid in pipes. To describe flow of fractal fluid between parallel planes, we should 
define new differential operators and the Navier-Stokes equations (1281) and (|3TT) . should be 
modified by using the product measure approach [63, [35J [64, [55] , 
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6 Fractal fluid with a r ^ 1 


Let us derive the Poiscuille equation for fractal fluids from the Navier-Stokes equations (l3Tj) . 
The Navier-Stokes equations for fractal fluid with a r = D — d 7 ^ 1 has the form 


.( 1 d 2 v x (r) d x +1 — a r dv x (r)\ 

A(d x ,a r ) ^ r 2 0r _ 2 g r 2 1 r 2a r -l J 


1 dp 
p dx 


0 , 


(48) 


where A(d x ,a r ) is defined by (1231) . d x = d — a x , and a x is dimension along the A"-axis. Using 
v x (r) as an effective scalar held p e ff(r) = v x (r), we can apply equations f[22|) . (12TT) and (1241) 
where D —$■ D x = D — a x and d —» d x = d — a x to get (148]) . Equation (1481) with a r = a x = 1 
gives ([35]) . 

For 1 < D < 3 and 0 < D < 1, the general solution of (HS|) is 


v x (r) = C\ r“ r dx + C 2 + 


dp 




2 (d x + a r ) a r A(d x , a r ) p dx 


(0 < D < 3, D^l). (49) 


For a r = a x = 1 equation (149|) gives (1361) . The constants C\ and C 2 in the general solution (1491) 
are determined by the boundary conditions 


v x (Ri) = v x (R 2 ) = 0. 


These conditions give the equations 
Ci R± r ~ dx + C 2 + 


d P O2a 


Ci R2 r ~ dx +C 2 + — - - — -— 5 ^ R z 2 a ‘ r = 0. 


2 ( d x + a r ) ay A(d x: a r ) p dx 
1 dp 


/?,, = 0 , 


->2 OL r 


Then the coefficients are 

Ci = 


2 (d x + a r ) a r A(d x , a r ) p dx 


dp R\ ar - Rl ar 


C,= 


2 (d x + a r ) a r A(d x , ay) p dx R 1 

1 


- R,?> 


dp R 2 ° r R° r ~ dx R^air ROlr—dx 


Substitution of 


2 (d x + ay) a r A(d x , ay) pdx R* r dx - R%’ 
and ([54]) into (149]) gives 

2(Xr T~)‘2iQ£r 


1 


v x (r = 


dp / R{ a '" - Rf 


2 (d x +a r )a r A(d x ,a r ) pdx\R* r dx - R% r dx 

r)2a r r>a r —d x p2a r z> 

it 2 J ^1 - 1 L 1 it' 


ty*<A-r d x j 


T)CX. r d x ryOc r d x 

LXY Ib 2 


‘lot.r T)Ot.r d x 

2 _ r 2a r 


(50) 

(51) 

(52) 

(53) 

(54) 


( 55 ) 


12 














for 0 < D < 3, where h^l. 

If R\ = 0 and f ?2 = R, then equation (T55l) has the form 

v (r) = _I_ ± R ** r ( (r_yr- d « _ /nH 

2 [d x + a r ) a r A(d x , a r ) fi dx \\R/ \R/ 

For a r — a x — 1 equation (15(1 gives (1551) . The cases a r < 1 and/or a x < 1 corresponds to 
fractal fluid. 

We can assume that a x > 1 can be used to describe fractal turbulent flow in pipe. This 
assumption is based on the fact that trajectories of the fluid particles are fractal curve, then 
a x > 1 (for example, the Koch curve with a x = ln(4)/ln(3) ~ 1.262). 



7 Fractal fluid discharge 


In general, fractal fluids cannot be considered as a fluid on fractal. Real fractal fluids have a 
characteristic smallest length scale such as the radius, Rq, of a particle (for example, an atom or 
molecule). In real fluids the fractal structure cannot be observed on all scales but only those for 
which R > Rq, where Ro is the characteristic scale of the particles. The concept of non-integer 
mass dimension of fractal fluid is based on the idea of how the mass of a fluid region scales with 
the region size, if we assume unchanged density. For many cases, we can write the asymptotic 
form for the relation between the mass Md(W) of a ball region W of fluid, and the radius R 
containing this mass as follows: 

M d (W) = Mo ( — J (57) 


for R/Rq 1. The constant M 0 depends on how the spheres of radius Ro are packed. The 
parameter D, which is interpreted as a dimension, does not depend on the shape of the region 
W, or on whether the packing of spheres of radius Ro is close packing, a random packing or 
a porous packing with a uniform distribution of holes. The non-integer mass dimension D of 
fractal fluid is a measure of how the fluid fills the integer n-dimensional Euclidean space it 
occupies. Note that the fact that a fluid is random or contains cavities does not necessarily 
imply that the fluid is fractal. 

Using the non-integer dimensional space approach, we can calculate the mass of fractal 
homogeneous fluids. Scaling law (157)) is obtained naturally in the framework of this approach. 
We can use the integration in a non-integer dimensional space [35] that is described by the 
equation 

drr D / dO sin D ~ 2 6 ip(r, 6), (58) 

Jo 

where d D r represent the volume element in the non-integer dimensional space. Using (158|) with 
<p(r, 9) — 1 , and 


IR D 


d u r<^(r) = 


27r (D-l)/2 re 

T((D - l)/ 2 ) J 0 


d6 sin D 2 6 = 


7T 1 / 2 r(U)/2- 1) 

r(D/2) 


(59) 
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( 60 ) 


we get the volume of D- dimensional ball with radius R in the form 

V t °/2 

Vd ^ T(D/ 2 + 1) RD ' 

The mass of fluid in W is described by the integral 

M d (W)= [ p(r) d D r, (61) 

Jw 

where r is dimensionless vector variable. For a ball with radius R and constant density p(r) = 
p = const, we get 

ty d / 2 a 

M d (W) = pV D = R ■ (62) 

This equation define the mass of the fractal homogeneous ball region of fluid with volume V/> 
For D — 3, equation (162)) gives the well-known equation for mass of non-fractal ball region 
M 3 = (4p7r/3) R 3 because T(3/2) = and T(z + 1) = zT(z). 

Let us determine the mass Q of fluid passing per unit time through any cross-section of 
the pipe (called the discharge). Not all pipe volume is occupied by fractal fluid. There are 
areas unoccupied by particles of fractal fluid. In continuum model of fractal fluids, we take into 
account this fact by using the integration in space with non-integer dimension 

2 7r d / 2 f R 

Q=p fml v ^ r)r ^ dr ' (63) 


where d = D — 1 is non-integer dimension of the cross-section, p is a constant density, and v x (r) 
is defined by equation ( 1461) . Substitution of (1461) into (1631) gives 


Q = 


p 7t( d l ^ 2 dp 

2 (D + 1) r((Z) — l)/2 ) pdx 


R d+1 


(0 < D < 3 D^l). 


(64) 


Note that p has physical dimension of mass (for example, kilogram). The mass of fractal fluid 
is thus proportional to (D + l)-power of the radius of the pipe. 

For D — 3, equation (1641) gives the well-known equation 


pvr dp 
8 p dx 


(65) 


The mass of non-fractal fluid is proportional to the fourth power of the radius of the pipe. 
The dependence of Q on dp/dx and R given by formula (l65|) was established empirically by G. 
Hagen in 1839 and J. L. M. Poiseuillc in 1840, and theoretically justified by G. G. Stokes in 
1845. 

Equation (164)) can be rewritten in the form 


PK d P r d +i 
8 Peff dx 


( 66 ) 
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with the effective dynamic viscosity 

Ml = ^ ^ (3_D,/2 r((D - l)/ 2 ) \i. (67) 

The mass (discharge) of fractal fluid is proportional to the non-integer power [D + 1) < 4 of 
the radius of the pipe, and the dynamic viscosity is effectively changed. 

For D — 3, equation (1671) gives /ieff = h- For 1 < D < 3, we have p, e ff > p. If 0 < D < 1, 
then fi e ff < 0. The effective dynamic viscosity of fractal fluid with D G (1; 3) increases 
with increasing a deviation of the dimension D from three. We can see an interesting effect 
of a negative effective dynamic viscosity for fractal fluid with dimension D G (0; 1). The 
strong fractality of the fluid, which is caused by small dimension, leads to an increased fluid 
flow compared with conventional medium. This is probably due to an increase in freedom of 
particles motion for the fractal fluid similar to 3D Cantor dust. 

The discharge function Q = Q(R, D) defined by (1641) for the different values of dimensions 
0 < D < 3 and the range R G [0,1] are present on Figures 5-8, where p = 1, p = 1, and 
dp/dx = —1. 

8 Conclusion 

In this paper, we propose a generalization of the Navier-Stokes equations to describe frac¬ 
tal fluids in the framework of continuum models with non-integer dimensional spaces. These 
equations contain generalized differential vector operators for non-integer dimensional space 
[72]173]. As an example of application of the suggested Navier-Stokes equations for fractal 
fluids, we consider a Poiseuille flow of an incompressible viscous fractal fluid in the pipe. The 
solution for steady flow of fractal fluid in a pipe and corresponding fractal fluid discharge have 
been derived. 

In this paper fractal fluid is described as a continuum in non-integer dimensional space. 
We assume that suggested continuum models with non-integer dimensional spaces and the 
correspondent Navier-Stokes equations for fractal fluids may be important for fractal theory of 
different type of media. 

As the main object for application of the proposed continuum models is a two-component 
medium, where distribution of one component (gas, liquid, solid) into another component (fluid, 
gas or empty space) can be characterized by non-integer mass or ’’particle” dimension. This 
non-integer dimensional component can be considered as a fractal fluid. One of the possible 
experimental methods for determining the presence of fractal properties of the two-component 
medium may be to use labels with radioactive isotopes for particles of component that is 
assumed a fractal. 

A basic idealized model of fractal fluid is a liquid distributed in empty space M 3 with non¬ 
integer mass dimension D < 3. In some sense the fractal fluid is considered as a liquid analog of 
fractal porous solid material. Fractal fluid can also be viewed as a two-phase medium consisting 
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of a liquid and a discharged gas instead of empty space, where the liquid is characterized by 
fractal mass dimension. 

As an object of study, we also can consider an emulsion, when both the dispersed and the 
continuous phase are liquids, and the dispersed phase is fractally distributed in continuous 
phase. An emulsion that is a mixture of two immiscible liquids, one of which (the dispersed 
phase) is fractally dispersed in the other (the continuous phase). In this case the dispersed 
phase can be described as a fractal fluid by suggested continuum models with non-integer 
dimensional space. The proposed models can be used for a solution that is a homogeneous 
mixture composed of only one (liquid) phase, where one phase has a fractal dimension. We 
can consider a fractal distribution of a solute dissolved in a non-fractal solvent, then the solute 
is considered as a fractal fluid. The solvent that is fractally homogeneously mixed with solute 
can be considered as a fractal homogeneous fluid. The homogeneity property of the fractal 
fluid means that two regions W\ and W 2 with the equal volumes V n {W\) = VniWz) have equal 
number of particles Nr,{Wi) = NniWz)- In other words the fractal fluid is called homogeneous 
if the power law Nr,{W) ~ R D (or Mo(W) ~ R D ) does not depend on the translation of the 
region W. 

We can consider a fractal distribution of small solid particles in the suspension. In this case, 
we have an internal phase (solid) that is fractally distributed through the external phase (fluid) 
by mechanical agitation. 

An object of investigations can be a liquid mixed with a solid particles, where the distri¬ 
bution of these particles in space can be characterized by non-integer mass dimensions, which 
can be caused by a power law distribution of particles by size or mass. 

As a complex medium which may exhibit fractal properties can be considered the blood 
that is composed of proteins, glucose, mineral ions, hormones, carbon dioxide, blood cells and 
other particles suspended in water. We assume that the blood as a multi-phase medium can 
have attributes of a fractal distribution for some blood components including bacteria, viruses 
and medicinal substances getting into the blood. 

We assume that the suggested approach to describe fractal fluids by continuum models with 
non-integer dimensional spaces may be important for fractal theory of blood flow in cardiovas¬ 
cular system, dynamics of fractal media in hydrologic modeling [66], [67], i68j [69] and it allows to 
develop the fractal dynamics of multi-phase media PHUZB- 
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Figure 1: Plot of the velocity function z = v(x,y) defined by (j45l) for the ranges x G [1; 100], 
y G [1; 100], and D — 2.9. 



Figure 2: Plot of the velocity function z = v(x,y ) defined by (|45l) for the ranges x G [1; 100], 
y G [1; 100], and D — 2.7. 
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Figure 3: Plot of the velocity function z = v(x,y) defined by (j45l) for the ranges x G [1; 100], 
y G [1; 100], and D = 2.0. 



Figure 4: Plot of the velocity function z = v(x,y ) defined by (|45l) for the ranges x G [1; 100], 
y G [1; 100], and D = 1.1. 
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Figure 5: Plot of the discharge function Q = Q(R, D ) defined by (|64l) for the range R G [0,1] 
and D G [2, 3]. 



Figure 6: Plot of the discharge function Q = Q(R, D ) defined by (EH) for the range R G [0,1] 
and D G [1, 2], 
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Figure 7: Plot of the discharge function Q = Q(R,D) defined by (j64li for the range R G [0,1] 
and D G [0.5,1.5]. 



Figure 8: Plot of the discharge function Q = Q(R, D ) defined by (1641) for the range R G [0,1] 
and D G [0, 3]. 
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